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In this letter we introduce the concept of stabilized vector solitons as nonlinear waves constructed 
by addition of mutually incoherent Townes solitons that are stabilized under the effect of a periodic 
modulation of the nonlinearity. We analyze the stability of this new kind of structures and describe 
their behavior and formation in Manakov-like interactions. Potential applications of our results in 
Bose-Einstein condensation and nonlinear optics are also discussed. 
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Since the introduction of the concept of soliton as soli- 
tary water waves with robust asymptotic behavior af- 
ter mutual collisions, many other physical systems have 
been found with similar dynamics, always described by 
nonlinear wave equations jlj. For solitons of Nonlinear 
Schrodinger equations (NLSE), the main interest in the 
early investigations was related with practical applica- 
tions in optical telecommunications, nowadays well es- 
tablished 1^. The recent interest on solitons in the field 
of Bose-Einstein Condensation (BEC) in alkali gases with 
negative scattering length |j, IH 1^1 1 shows the timeli- 
ness of the topic and its central place in modern physics. 

Despite the success of the concept of soliton, these 
structures arise mostly in l-|-l-dimensional configura- 
tions. In the NLSE case this is mainly due to the well 
known collapse property in multi-dimensional scenarios 
0. In the Optical context, collapse means that a laser 
beam with power higher than a critique threshold, will be 
strongly self-focused to a singularity when propagates in 
a Kerr-type nonlinear medium, whereas for lower powers 
it will spread as it propagates. This behavior has also 
been observed in experiments with matter waves . 

Since collapse prevents the stability of multidimen- 
sional "soliton bullets" in systems ruled by NLSE, a great 
effort has been devoted to search for systems with sta- 
ble solitary waves in multidimensional configurations . 
A new way to generate stabilized two-dimensional soli- 
tary waves has been recently proposed for optical systems 
|lCllll| . The idea is to prevent collapse by using a spatial 
modulation of the Kerr coefficient (the nonlinearity) of 
the optical material so that the beam becomes collapsing 
and expanding in alternating regions and is stabilized in 
average. The idea has been extended to the field of mat- 
ter waves in Refs. 0|. Finally, in Ref. |Q some 
general results are provided. 

In the present paper we will extend this analysis to the 
case of mutually incoherent beams with unexpected and 
surprising results. This is, to our knowledge, the first 



theoretical evidence of two-dimensional stabilized vector 
solitons (SVS), a new kind of nonlinear waves which can 
be constructed in two ways: by direct combination of 
several Townes solitons or as a result of Manakov inter- 
actions [III between Townes solitons. In both cases, the 
stabilization against collapse is obtained by the effect of 
a peridodic modulation of the nonlinearity. 

The model.- Let us consider a n-component system 
modelled by equations of the type 




(1) 



where j = 1, . . 

I dy^ ,ajk G IK are the nonlinear coupling coefficients 
and g(t) is a periodic function accounting for the modu- 
lation of the nonlinearity. 

E qs. l H) are the natural extension of the Manakov sys- 
tem jl5| to two transverse dimensions and an arbitrary 
number of components. In Optics, for spatial solitons t 
is the propagation coordinate and Uj are n mutually in- 
coherent beams. One-dimensional Manakov-type models 
have been extensively studied in nonlinear optics, mainly 
due to the potential applications of Manakov solitons in 
the design of all-optical computing devices ^16J. In BEC 
these equations (with an additional trapping term) de- 
scribe the dynamics of multicomponent condensates, Uj 
being the wavefunctions for each of the atomic species 
involved [Ulii. 

In the scalar case (n = 1), it is well known that, 
if g is constant, there is a stationary radially symmet- 
ric solution of Eq. ^ (the so-called Townes soliton): 
M(r, t) = $(r)e''*'*. This solution is unstable meaning that 
the value of the norm N^^^ = {J |<I>p)-^/^ is critical in the 
sense that any generic slight perturbation of the initial 
condition will yield to collapse or spreading of the dis- 
tribution. We must also notice that, due to the scaling 
invariance of the cubic NLS, a family of Townes solitons 
can be generated by making ^\{r, t) = A<I>(Ar, X^t). 
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FIG. 1: [Color online] (a) Stabilized vector solitons for 
Qfi = Qf2 = l/\/2 (the evolution for U2, not shown here, is 
very similar). Shown are the time evolution of the width 
Wi = (/|{2 (s^^ + J/^)!^'!^)"'"''^ and amplitude max^gg2 \ui{x,t)\. 
The insets show pseudocolor plots of \ui{t, x, y)\^ for different 
times, (b) Same as (a) but for ai — l/\/3, Q!2 = ^2/3. For 
both cases g{t) = —2n + 87rcos(40t). 



It has been shown in 'l4f| that an adequate modulation 
of the nonlinearity, will asymptotically stabilize a Townes 
soliton yielding to a rapidly oscillating stabilized Townes 
soliton (STS), which we refer hereatfer as ^s- In this 
paper, we take g{t) = go + gi cos fit but we expect that 
most of our results with similar periodic functions will be 
qualitatively the same 14] . 

Stabilized Vector Solitons (SVS).- For a given set of 
parameters ajk it is possible to use stabilized Townes 
solitons to build explicit solutions of Eqs. These so- 
lutions are constructed by taking Uj = $5 = Q;j$5,j = 
1, 



, n for any set of coefficients aj satisfying 
ajitti + ... + ajna^n = 1; j = 1, • • ■ , 



(2) 



It is not obvious that these new solutions will be sta- 
bilized by a periodic modulation g{t) of the nonlinear- 
ity. Writing Uj = $5^- -|- 5j the equations for 5j contain 
cross-modulation terms which could lead to growth of 
these small perturbations. To test that the stabilization 
is possible in a wide range of configurations, we have 
considered several important examples. First, we have 
studied the most relevant case n = 2 and integrated nu- 
merically Eqs. 1^ with different initial data of the form 
Uj = aj^s satisfying |(5J) and found that these new vec- 
tor solitons remain stabilized as shown in Fig. ^ From 
now on we will name these structures as Stabilized Vec- 
tor Solitons (SVS). We will see below how they emerge in 
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[Color online] Fast collisions of stabilized Townes 



FIG. 2: 

solitons. Initial data are iii(0,r) — e'^i''<l>(|r + rij), 
U2{0,r) = e'"2--$(ir-Hr2]) with vi = (5/^2, 5/v^), va = 
(-5/a/2,5/^) and n = (-6,-6), ra = (6,-6). (a) Sur- 
face plots of \ui\'^ and |u2p for different times, (b) Evolu- 
tion of the widths VFi,^(t) = (/ (a;- < x >f \ui{x,y,t)\^Y'^, 
Wi,y{t) = {J{y- < y >f\ui(x,y,t)fy/^. (c) Evolution of 
the maximum amplitude maxj^, y)gu2 \ui{x,y,t)\. 



collisions of mutually incoherent STSs, which correspond 
to the so-called Manakov interactions 15]. 

We have studied other situations such as a symmetric 
superposition of four STS (n — 4) with aj = l/V^ and 
found similar results. Thus, these structures exist in a 
wide range of parameters and configurations. 

Manakov interactions of STS.- Depending on the mu- 
tual velocity of the two interacting STS, we have divided 
the regime of collisions in two different ones: fast and 
slow collisions. As we will see below, one of the main re- 
sults of our work is the possibility of obtaining SVS after 
slow collisions of STS. 

First we have studied collisions of "fast" STSs after 
which the solitons emerge with only moderate modifica- 
tions of their amplitude and width as is shown in Fig. 
El It can be seen [Fig. Efb)] that during the collision, 
the soliton becomes spatially asymmetric. An internal 
asymmetric breathing mode of small amplitude is excited 
which decays at longer times (not shown in the figure) to 
the "normal" symmetric breathing mode shown by STSs. 

These behaviors can be accounted for by a finite- 
dimensional reduction of Eqs. ^ by means of the time- 
dependent variational approach. Notice that Eqs. ^ 
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can be obtained from the Lagrangian density 



I ( dul 



dt 



(an 



dt 

litil"^ + 2ai2|uinu2p + a22|u2|'^) . 



We choose a simple ansatz accounting for head-on sym- 
metric coUisions of equal stabilized solitons moving with 
opposite speeds and centered on (— i?, 0) and 0) 

Although gaussians do not have the right asymptotic de- 
cay as STSs, our choice simplifies the calculations and is 
enough for our present objectives. The standard varia- 
tional method [l9| leads to the equations (for ajk ~ 1) 
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together with the complementary relations Pj — Wj /2Ljj^ 
{j — x,y), V — £~2£/3x, and the conservation law N{t) — 
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TT\A(0)\^u}xC^)oJy(0). The different terms in 



Eqs. I@J account for the phenomenology shown in Fig. 
Inland other "fast collisions" studied. For example, they 
contain the asymmetric interaction (notice the differences 
between Eqs. <j4b|l and lBc|l ') due to the fact that both 
solitons approach along the x axis and thus become more 
elongated along that direction as seen in Fig. Hfb). We 
have numerically integrated Eqs. Q for fast collisions 
taking as initial data stabilized gaussian functions [lllll^ 
and have found results similar to those shown in Fig. [3 
The regime of slow collisions is in the range |t'2 — t^i | ^ 
3. In this case the collisions of STSs lead to the for- 
mation of two vector solitons as shown in Fig. |SIa,b). 
It is remarkable and one of the main results of the pa- 
per that the collision mechanism allows the complex co- 
herent rearrangement necessary for the formation of the 
vector solitons. The fraction of "mass" interchanged 
by the incoming solitons is a function of the only rel- 
evant parameter for direct collisions \v2 — vi\ (due to 
the Galilean invariance) as shown in Fig. EJa). In 
the range 0.2 < 11)2 — I'll < 3 we observe formation of 
two vector solitons which seem to be either unstable or 
performing high amplitude oscillations for higher speeds 
and stable in the lower range of speeds (approximately 
0.2 < |w2 — wi| < 1.2). If the speed is decreased further we 
observe two outgoing vector solitons with complex tran- 
sient dynamics and nontrivial dependence of the fraction 
transferred as a function of 1112 — 




FIG. 3; [Color online] Head-on collisions of STSs for initial 
data v\ — —V2 = 0.3. (a) Surface plots of \u\\^ , (b) sur- 
face plots of \u2\^ ■ (c) Evolution of the maximum amplitude 
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FIG. 4: Asymptotic behavior after head-on collisions of STSs 
of the form m = $(|r - £|)e™i^, Ma = "I-dr + ^|)e*"2== and 
large enough ^ (~ 4). (a) Regimes of behavior as functions 
of vi,V2. (b) Quotient (/) of the squared amplitudes of the 
small and large peaks which are generated after the collision 
when a VS is formed [see Fig. |2] for the regime of speeds in 
the range 0.2 < |«2 - «i| < 1.2. 



Finally, if the initial speed of the colliding solitons is 
very small or zero, we have observed a quasi-bound state 
of two SVSs which shows several recurrent collisions as 
shown in Fig. |31 From our simulations we cannot con- 
clude if this is a true bound state or it finally decays to 
vector solitons. In Fig. ^b) we summarize the results of 
our numerical exploration of STS collisions. 

We want to point out that Eqs. Q provide a reason- 
ably good description of the phenomena described here 
as far as the ansatz given by Eqs. Q can describe these 
complex dynamical behaviors. An example: for very low 
speeds the variational equations predict the formation of 
an oscillating bound state of two STSs. Although this 
is not the real behavior (a bound state of two SVSs is 
formed) we get bound states. 

The formation of vector solitons from stabilized scalar 
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FIG. 5: [Color online] Oscillations after a collision of STSs 
with vi = V2 = 0. (a) Isosurface plot of \ui{t,x,y)\'^ 
(shownare isosurfaces corresponding to 0.05, 0.15 and 0.25), 
(b,c,d,e) Surface plots of \ui{t, x,y)\^ for (b) t = 0, (c) t = 50, 
(d) t = 100, (e) t — 150. The corresponding evolution for U2 
is symmetric with respect to the y axis. 

solitons is a nontrivial phenomenon since there is a del- 
icate balance of both components which must be satis- 
fied in order to avoid destabilization either to collapse 
or expansion of these structures. It is curious that the 
system is able to interchange just the right amount of 
energy to keep both solitons bounded. In fact, the colli- 
sion mechanism described here can be seen as a way to 
generate appropriate stabilized vector solitons up from 
STSs which could be otherwise difficult to obtain. This 
is another proof of the structural stability of these new 
physical objects. 

It is also remarkable that no collapse phenomena is 
observed in our simulations, rather instead most of the 
collisions observed lead to remarkably robust scalar or 
vector solitons. This is very different from what happens 
in coherent collisions of STSs which lead to collapse. 

Conclusions.- In this paper we have described a new 
type of vector solitons, the stabilized vector solitons 
(SVS). We have studied their stability and shown how 
they arise in Manakov collisions of STSs. Other phe- 
nomena seen in collisions of STSs have been described 
and analyzed such as "quasi" -bound state formation and 
excitation of asymmetric oscillation modes. 

We must stress that the limit n — > oo of our model 
could be used to study nonlinear propagation of totally 
incoherent light. In BEG a detailed investigation of the 
above system could have implications in the study of de- 
coherence effects in BECs. 
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